Strong magnetoresistance induced by long-range disorder 
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We calculate the semiclassical magnetoresistivity pxx{B) of non- interacting fermions in two di- 
mensions moving in a weak and smoothly varying random potential or random magnetic field. We 
demonstrate that in a broad range of magnetic fields the non-Markovian character of the transport 
leads to a strong positive magnetoresistance. The effect is especially pronounced in the case of a 
random magnetic field where pxx{B) becomes parametrically much larger than its _B = value. 
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The magnetoresistance (MR) is one of the most fre- 
quently studied characteristics of the two-dimensional 
electron gas (2DEG). When the effect of disorder is de- 
scribed by a collision integral within the semiclassical 
Boltzmann equation approach, the resistivity tensor p{B) 
for an isotropic system has the Drude form 



m = 



^-1 



(1) 



where n is the carrier density, m the effective mass, 
UJc = eB/mc the cyclotron frequency, and r the trans- 
port scattering time. In particular, the longitudinal 
resistivity pxx is independent of the magnetic field B, 
Pxx{B) — po = m/e^riT, irrespective of the form of the 
impurity collision integral. This result is solely deter- 
mined by the Markovian character of the transport as- 
sumed in the Boltzmann equation description. 

Deviations from the constant Pxx{B) are convention- 
ally termed a positive/negative MR, depending on the 
sign of the deviation. The negative MR |0] induced by the 
suppression of the quantum interference correction by the 
magnetic field is a famous manifestation of weak local- 
ization. Another source of negative MR is the Altshuler- 
Aronov correction to the conductivity due to enhance- 
ment of the electron-electron interaction by the diffusive 
motion of particles 0]. Both these effects are of quan- 
tum nature and lead to a correction of order /h to the 
conductivity <Txx, and thus to a small correction to pxx- 

However, as we will show, already at the classical level 
there exists a non-trivial MR which can be much stronger 
than the quantum one, if the correlation length d of disor- 
der is sufficiently large, kpd ^ 1 (where kp \s the Fermi 
wave vector). This is due to memory effects which are 
neglected in the collision integral description of disorder. 

Transport properties of the 2DEG in a smooth ran- 
dom potential (RP) V{y) are of particular interest, since 
in currently fabricated high-mobility semiconductor het- 
erostructures the disorder has long-range character. The 
high mobility of these samples is achieved by placing the 
charged donor ions in a layer separated by a large dis- 
tance d {kpd 10) from the 2DEG plane. Assuming 



the positions of these impurities to be statistically dis- 
tributed with a sheet density rii, the correlation function 
WV(r — r') = {V{r)V{r')) is given in momentum space 

by 



-2qd 



(2) 



A new type of transport problem occurs in these sys- 
tems when a large magnetic field B 2± i?i/2 — 2{hc/e)n 
is applied such that the lowest Landau level is approxi- 
mately half filled. The metallic state then observed has 
been described |^ in terms of composite fermions (CF's) 
moving in a weak effective field B = B — Bi/2- The 
CF's are scattered by an impurity-induced random mag- 
netic field (RMF) 6B{r) characterized by the correlation 
function WB(r — r') = {5B{r)5B{r')) with Fourier com- 
ponents 



(3) 



While the above RMF is fictitious, a real long-range cor- 
related RMF can also be realized in semiconductor het- 
erostructures by attaching superconducting or fer- 
romagnetic overlayers. We will study below the 
case of a weak RP or RMF, which means that I ^ d, 
where I is the mean free path in zero average field B — Q. 
Let us stress that we consider a situation with only the 
smooth disorder (RP or RMF) present. This should be 
contrasted with the starting point of Q , where the resis- 
tivity was assumed to be dominated by a white-noise RP 
while a weak long-range RMF was considered as a small 
perturbation. 

As one manifestation of the strongly non-Markovian 
transport, it has been shown recently |]8Hl0| that in a 
sufficiently strong B the MR drops exponentially with 
B because of a "classical localization" effect caused by 
adiabacity of the motion. This holds true both for the 
motion in a RP g and in a RMF The condition 

of this adiabatic regime is ujcj ^ (l/dy^^. 

In this paper we study the region of smaller magnetic 
fields in which different non-Markovian processes become 
important. We find that the exponential falloff of pxx is 
preceded by a strong positive MR. The effect is especially 
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pronounced in the case of the RMF, where the increase 
of pxx is much larger (in the weak disorder Umit) than 
its zevo-B value po = m/e^riT. 

We first outline the physics of the effect on a qualita- 
tive level. The zeio-B transport scattering rates in the 
RP and RMF are given by Q 



dqq^Wviq) (RP); (4) 
dqWeiq) (RMF). (5) 



e \2 1 
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Let us now discuss the nature of the particle trajectories. 
The change d(j> of the polar angle of the particle velocity 
in the time interval dt is given by 



d(P : 



{mvFy^n^e,jdjV{r)dt (RP); 
{e/mc)SB{v)dt (RMF), 



(6) 
(7) 



where n = v/|v| is the unit vector in the direction of the 
velocity. Taking into account that the trajectory within 
the correlation domain of the size ~ d is almost a straight 
line, one gets (^, (^). One source of the MR is the bend- 
ing of the trajectory within the correlation domain by 
the magnetic field B. This leads to a small negative 
MR, of the relative magnitude - {d/Rc)^ < 1 
There exists, however, a much stronger effect related to 
returns of the particle to spatial regions close to the start- 
ing point. Under the condition lUcT ^ 1 the particle 
trajectory is a sequence of slightly distorted cyclotron 
circles. The center of the orbit is shifted by a random 
vector 6 after one cyclotron revolution, with {S) = 0, 
{SD = (Sy) = 2ttP /{lUct)^ ■ If, for definiteness, the veloc- 
ity v||x at the time t — 0, the trajectory passes a distance 
6y from the original point at the time t ~ 2'k/uJc- The 
correlation of the scattering processes (^), (0) at i ~ 
and t ~ 2n jiOc induced by the long-range character of 
the disorder leads to a correction to the relaxation rate 
and, therefore, the resistivity 



[d/{8Y'^f (RP); 

rf/(^2^1/2 (RMF). 



(8) 
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The same is valid for t ~ ^T^njuc (the n-th cyclotron 
revolution, n = 2, 3, . . .), with (J^) multiplied by n. For 
the RP case the corresponding sum over n converges, 
^^-^ = C(3/2), leading simply to the renormal- 

ization of a numerical factor in (^. In fact, in this case 
the above consideration can be made fully quantitative, 
and the result agrees with what we will find below from 
the Liouville equation (|lO|). The corresponding positive 
MR, Ap-cx/po ^ {d/lY{ijJcTY^'^ ^ becomes of order unity 
at the upper bound of the considered range of the mag- 
netic fields, uicT < (l/d)"^^^. Here the system enters the 
adiabatic regime. 

In the RMF case the sum over n is of the form ^ n^^/^ 
and is thus determined by the upper cutoff, which is n ~ 



ujcT. The resulting correction Ap^x/po {d/l){(-OcTY 
reaches a value of order unity at lOcT ~ (Z/d)^/^, i.e. 
far from the adiabatic regime. At larger magnetic fields, 
where Ap^x/ Po ^ 1, a self-consistent treatment is needed 
(see below). 

We describe now a formalism which allows us to calcu- 
late the MR more systematically. We consider first the 
RMF case. The starting point is the Liouville equation 



(io + 5L)g{w,Y 



cos ( 



mc 0(p 



SL 



^e)\ (10) 



for the deviation Sf{uj,r,(j>) = eEvF^^g{oJ,r, cj)) from 
the equilibrium distribution function /o — 0{eF — e). 
Here E = i? (cos (/<£;, sin (^£;) is the electric field and 
n — (cos (j), sin 0) the unit vector determining the ve- 
locity direction. The current density is given by j = 
~e J (^j^^vJ/, which yields the conductivity tensor 




(Lo + SL)- 



COS ( 

siuf 



(11) 



where Np is the density of states and the angular brack- 
ets denote the averaging over configurations of the RMF 
SB{r) with the correlation function (^). 

Expanding ( p] ) in SL, averaging over the disorder and 
resumming the series (in the same way as it is done for 
a quantum-mechanical Green's function), one arrives at 



{Lo + Af); 



(12) 



Ln = 



— VjJ OJc 

—LOc —iuj 



where M is the "self-energy" and the 2x2 matrix M is 
defined by 



M 



cos ( 
sin f 



M 



COS( 

sine 



(13) 



In zero B the self-energy M is given in the leading 
approximation by the first term of the perturbative ex- 
pansion, M = —{5L Lq"^ 5L), yielding 



Mrr = -2« 



X smt 



6x2/" d'^q 



mc) J (27r)2 2tt 

vpq cos {4> — 4>q) — uj 



smc 



(14) 



where (j)q is the polar angle of the momentum q. Taking 
into account that u should have an infinitesimal posi- 
tive imaginary part [uj lo + iO) and considering the 
limit w ^ 0, we get M^x = 1/t with 1/r given by (j|). 
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Equation (|1J) reproduces then the Drude resistivity, as 
expected. 

Now we calculate the B-dependent correction to jl^ ) 
determined by the return processes described above. For 
this purpose, we have to replace the free propagator 
Lq^ = {—iuj + ivpq cos {(j) — (f>q))^^ entering ( p^ by the 
one describing the motion in the magnetic field in the 
presence of disorder. At small q this would be simply the 
diffusion propagator; however, since we are interested in 
the short-scale physics, q > Rc^i the diffusion approxi- 
mation is not appropiate. Using the fact that the particle 
is only scattered by a small angle within the correlation 
length d, we can approximate the motion by a Fokker- 
Planck equation corresponding to the diffusion in mo- 
mentum space: 



AM,., = 2(— )' 
■.mc/ 



(Pq 



(2^)2 

X / ^ sm0gD(w,q, 



Weiq) 



(15) 



where gj^ is the solution of 

—iLO + ivq COS 



d 1 52 

' d(f) T d(j) '^ 



5D(t^,q, 0) 
(16) 



Now we make use of the fact that the solution of jl^ ) 
determines the spatial dispersion of the conductivity in 
the situation of small-angle scattering : 

ayy{uj,q) = e^vlNp — sin gD(w, q, </>) . (17) 

At a; = one has cryy(a\\y) = For q||x we get, 

using the solution of ( p^ ) at lUcT ;» 1 JTif , 



2ao r^d<P 



LOcT Jq 2tT 



y sin (A sin (/-'e^-^^'^''^') 



\G^Jl{qR,)l{qR^f , 



(18) 



where K[(^, 0') — iqRc{sm(j) — sin cf)') — [{qRcY /2LJcT]{(j) — 
(j)') and (To = Pq^- This determines the integrand of jl^ ) 
for general orientation of q, with the result for lo — 



Kmc) ttW^ Jq q 



dq j2 



JtiqRc)WB{q) 



/ e \2 T 



Wb{0) 



(19) 



In agreement with the above qualitative picture, the main 
contribution to the integral comes from q R^^ ■ There- 
fore, we neglected the g-dependence of Wsiq) in the sec- 
ond line of (^, in view of d < i?c. The MR is thus 
positive and quadratic in B, 

^P..Ip^ = {B/Bof = 4a^u,T)^ = 2{d/l)(uj,Tf , (20) 



where Bq = {5B'^{r)}^/^ = Wg^^{0) is the amplitude 
of the RMF fluctuations and a ^ 1 is the parameter 



characterizing the RMF strength ||,|l^, a — dcuc / 



VF 



with 



,(0) 



= eBa/mc. Eq. (pd) is valid for B Bq, 



whereas the adiabatic regime begins at B ~ B^a 
In the intermediate range, Bo < B < the 
positive MR gets large, IS.pxxl Pa 3> 1. In this region, 
Eq. ( p9| ) should be treated self-consistently, i.e. r in the 
r.h.s. should be understood as a renormalized scattering 
time, = Tg"^ + l^M^x- The result for B > is 
Pxx/pa = B/Bq, or, in the form valid both below and 
above Bq, 



Pxx/po = l/2+[l/A+{B/Boy] 



211/2 



(21) 



which is our main result for the case of RMF. At 
B/Bq a^^/^ the resistivity reaches its maximum 
Pxx ^ ct^^^^Po- in still higher fields pxx drops rapidly 
due to the adiabatic character of motion. 

We turn now to the RP scattering. The operator SL 
in (not) has then the form 



d 

5L = 6v{r)ny + {y6v{r))n± — , 



(22) 



where Sv{r) — V{r)/pF is the spatial variation of the 
Fermi velocity and nj^ = z x n = (— sin 0, cos 0). At zero 
B we get, instead of (|lj). 



Pf 



d^q 



f27r)2 2n 



fWv{q) 



sm sm 



.) 



■ sm sm ( 



</-?), (23) 

VFqCOS[0 ~ 0qj — LU 

reproducing the result (^. To calculate the MR, we re- 
place, as in the RMF case, the free propagator in ( [2^ ) 
by ijj^, where Ld is the operator in the square brack- 
ets in (p^. Solving again the equation Lugu{oJ,(l,(j>) = 
sin (/) sin (0 — 0,) at ^ 1, we find after some algebraic 
manipulations 



AM, 



1 



2'KPpVF 



dqq^Wviq) 



coth 



TTjqRc)' 

2uJcT 
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In contrast to (p^), the integral is determined by the re- 
gion of momenta q^ ^ uJct/R^ ^ (S^)^-^, as might be 
expected from the above qualitative consideration. Since 
d <^ {S^y^^, we can again neglect the g-dependence of 
Wv{q), which yields for UcT < {l/d^/^ 



Apxx _ C(3/2) rWyjO) , 
po ' 47r2 p|,t;^Z3 

2C(3/2) (d^^ ' 



(uJct) 



9/2 



9/2 



(24) 



We have performed numerical simulations of the MR 
for both types of disorder (RMF and RP) . In Fig. |l| the re- 
sults for the RMF are shown, for three different strengths 
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of the disorder (a = 0.2, 0.083, 0.0138). At a < 1, 
the theoretical prediction of the strong positive MR ( |2l| ) 
crossing over to the negative one at _B '-^ Boa~^^"^ is fully 
confirmed by the data. At moderately small a {a = 0.2 
in Fig. |l|) the positive MR still exists, but becomes weak; 
this is the region of a relevant to the composite-fermion 
description of the vicinity of half filling of the lowest Lan- 
dau level {i/ — 1/2). The numerically calculated MR for 
a ^ 0.2 — 0.3 agrees well 10 1 with the experimental data 
around i' = 1/2. At large a > 0.5 the region of posi- 
tive MR disappears, and pxx drops monotonously with 

B m. 
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FIG. 1. Magnetoresistivity [normalized by the Drude value 
po = m/e^riT with 1/r given by ^] in random magnetic field 
from the numerical simulations for three diff'erent strengths 
of the disorder; the full line corresponds to Eq. 




FIG. 2. Magnetoresistivity in a random potential from 
computer simulations in comparison with Eq. (124 



The numerically found MR for the RP case (Fig, p) 
shows good agreement with the theoretical result (24) 
up to Apxx/po 1- At larger B, p^x deviates from (24) 
and starts to decrease, as expected. 

In conclusion, we have demonstrated that the 2D 
fermion gas shows for weak long-range correlated disor- 
der a strong positive MR in moderately strong magnetic 
fields 1 <^ oJcT < (l/d)"^^^, due to the non-Mar kovian 
character of transport. The effect is especially pro- 
nounced in the case of the RMF. Our findings explain, 
in particular, the positive MR of composite fermions ob- 
served experimentally around u = 1/2. 
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